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Q ■ Abstract 

\o : 

We pursue the study of SU(2) Euclidean Yang-Mills theory in the maximal Abelian gauge 
by taking into account the effects of the Gribov horizon. The Gribov approximation, previ- 
ously introduced in pQ, is improved through the introduction of the horizon function, which 
Qh. is constructed under the requirements of localizability and renormalizability. By following 

Zwanziger's treatment of the horizon function in the Landau gauge, we prove that, when cast 
in local form, the horizon term of the maximal Abelian gauge leads to a quantized theory 
which enjoys multiplicative renormalizability, a feature which is established to all orders by 
means of the algebraic renormalization. Furthermore, it turns out that the horizon term is 



compatible with the local residual U(l) Ward identity, typical of the maximal Abelian gauge, 

1/2 

which is easily derived. As a consequence, the nonrenormalization theorem, Z g Z^ = 1, re- 
lating the renormalization factors of the gauge coupling constant Z g and of the diagonal 
gluon field Za, still holds in the presence of the Gribov horizon. Finally, we notice that 
a generalized dimension two gluon operator can be also introduced. It is BRST invariant 
on-shell, a property which ensures its multiplicative renormalizability. Its anomalous dimen- 
sion is not an independent parameter of the theory, being obtained from the renormalization 
factors of the gauge coupling constant and of the diagonal antighost field. 
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1 Introduction 



The maximal Abelian gauge j2j El Ej provides a useful tool in order to investigate nonperturbative 
aspects of Yang-Mills theories. This gauge turns out to be suitable for the study of the dual 
superconductivity mechanism for color confinement [5J, according to which Yang-Mills theories in 
the low energy region should be described by an effective Abelian theory [01 |HJ E] in the pres- 
ence of monopoles. The condensation of these magnetic charges leads to a dual Meissner effect 
resulting in quark confinement. Here, the Abelian configurations are identified with the diagonal 
components A 1 , % = 1, N — 1, of the gauge field corresponding to the (N — 1) generators of 
the Cartan subgroup of SU(N). The remaining off-diagonal components A a , a = 1, N 2 — N, 
corresponding to the (iV 2 — N) off-diagonal generators of SU(N), are expected to acquire a mass 
through a dynamical mechanism, thus decoupling at low energies. 

The maximal Abelian gauge is a renormalizable gauge in the continuum ^01 EU 112 Ej, while 
possessing a lattice formulation [Hill]. This property, shared only by few other gauges, such as the 
Landau and Coulomb 1 gauges, gives to the maximal Abelian gauge a privileged role. It allows for 
a comparison between the results obtained in the continuum and in lattice numerical simulations. 

Numerical studies of the behavior of the gluon propagator in the case of SU(2) can be found 
in In particular, according to both diagonal and off-diagonal gluon propagators turn 

out to be suppressed in the infrared region. The diagonal propagator exhibits a Gribov like be- 
havior, while the off-diagonal transverse component is of the Yukawa type. The value reported for 
the off-diagonal gluon mass is of approximately 1.2GeV fH>\ HE], and turns out to be twice bigger 
than the mass parameter entering the diagonal propagator. As a consequence, the off-diagonal 
propagator is short-ranged as compared to the diagonal one, a feature which is in agreement with 
the Abelian dominance hypothesis [01 El E] • 

As the Landau and Coulomb gauges, the maximal Abelian gauge is affected by the existence 
of the Gribov copies ^Zj 2 , which have to be taken into account in order to properly quantize the 
theory. A first step in this direction was achieved in pQ, where the original construction outlined 
by Gribov for the Landau gauge ^7j was generalized to the maximal Abelian gauge. In particular, 
we have been able to generalize to the maximal Abelian gauge Gribov's result stating that for any 
field close to a horizon there is a gauge copy, close to the same horizon, located on the other side 
of the horizon ^7j. This result has provided support for restricting the domain of integration in 
the Feynman path integral to the so-called Gribov region, i.e. to the region in field space whose 
boundary is the first Gribov horizon, where the first vanishing eigenvalue of the Faddeev-Popov 
operator appears. 

According to Gribov's original procedure [T7j, this restriction can be implemented by means of a 
no-pole condition on the ghost two-point function, which can be worked out order by order [T7| . 
In particular, the so called Gribov approximation amounts to work out the no-pole condition at 
the first nontrivial order, resulting in a new nonlocal term which has to be added to the Yang- Mills 

^dn the case of the Coulomb gauge we should remind that the issue of the renormalizability is still under debate, 
see 23] for a detailed account. 

2 See 2H| f° r the construction of an explicit example of a zero mode of the Faddeev-Popov operator. 
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action. In the Landau gauge, this term reads [T7j 

S G = -g 2 i A N J d*xA A ±A A , (1) 

where the index A belongs to the adjoint representation of SU(N), A — 1, N 2 — 1, and where 
7 is the Gribov parameter [Uj. It is not a free parameter of the theory, being defined by a gap 
equation which enables us to express it in terms of the gauge coupling g. At the first order, the 
gap equation defining 7 is given by 

Although expression (JTJ), being quadratic in the gauge fields A A , is useful to analyse the modifica- 
tions of the tree level gluon propagator due to the restriction to the Gribov region, it suffers from 
several limitations. For instance, it does not allow to go beyond the tree level approximation, by 
taking into account in a consistent way the effects of the higher order quantum corrections, which 
would provide a better comparison with the results obtained from lattice numerical simulations 
as well as from the studies of the Schwinger-Dyson equations. The resolution of this important 
issue is due to Zwanziger who, by a recursive characterization of the lowest eigenvalue of the 
Faddeev-Popov operator, —d^D AB = — <9 M (d^S AB + gf ACB A^) , has obtained a closed expression 
for the nonlocal term which implements the restriction to the Gribov region in the Landau gauge 
[TTH I2U] . This term, known as the horizon function, is given by 

Szw = -</V / *x f ABC A B [{d,D^] AD f DEC A E . (3) 

Notice that Zwanziger's horizon term, eq.Q, reduces to the Gribov term {IJ in the quadratic 
approximation, in which the Faddeev-Popov operator, —d^D AB , is replaced by the Laplacian, 
i.e. by —S AB d 2 . To some extent, Zwanziger's solution amounts to perform a re-summation of 
all higher order terms neglected in Gribov's quadratic approximation. In spite of its apparent 
nonlocality, the horizon function, eq.®, can be cast in local form through the introduction of 
a suitable set of additional fields. Remarkably, the resulting local action turns out to be multi- 
plicatively renormalizable to all orders ^01 1201 IS] • This property has made possible to evaluate 
the higher orders quantum corrections by taking into account the effects of the restriction to the 
Gribov region. In [22] one finds a study at one loop order of the gluon condensate (A A A A } in the 
presence of the horizon term Q. More recently, the evaluation of the two loop quantum correc- 
tions to the gluon and ghost propagators and their relationship with the infrared freezing of the 
running coupling constant has been worked out in [23 121] • We emphasize that the restriction to 
the Gribov region, as implemented by the Zwanziger horizon function, gives us a useful theoretical 
framework in order to achieve a rather clear understanding of the origin of the infrared suppression 
of the gluon propagator as well as of the infrared enhancement of the ghost propagator observed 
in the Landau gauge in both lattice numerical simulations [23 1211 123 EU 123 EEH EH E21 EH] and 
Schwinger-Dyson approach jHU ESI EH EH EH] • 

Concerning now the maximal Abelian gauge, till now, the restriction to the Gribov region has been 
implemented within Gribov's quadratic approximation, as reported in [I] in the case of SU(2). 
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For the analogue of the nonlocal term (JI} we have obtained 



S%™ = -2rf? JdhA^A^ (4) 
where is the diagonal component of the gauge field, defined according to the decomposition 

A, = A^T a + A,T 3 , (5) 

where T a , a = 1,2, denote the off-diagonal generators of 577(2), while T 3 stands for the diagonal 
generator, 

[T a ,T b ] = ie ab T 3 , 

[T 3 ,T a ] = ie ah T\ (6) 

where 

£ ab _ _ £ ab3 

£ ab £ cd = 6 ac 6 bd _ 6 ad 6 bc _ ^ 

As in the case of the Landau gauge, the Gribov parameter 7 is not free, being determined pQ by 
a gap equation 3 similar to eq.(J2J)- The addition of the Gribov term to the Yang-Mills action 
deeply modifies the infrared behavior of the gluon and ghost propagators in the maximal Abelian 
gauge. Let us mention here that the results which we have obtained |Tj for the diagonal and off- 
diagonal components of the gluon propagator are in qualitative agreement with those reported in 
[IS], see next section for a brief review. Nevertheless, the nonlocal term (J3J) suffers from the same 
limitations of the corresponding term, eq.(fTj). of the Landau gauge, requiring thus the construction 
of the analogue of Zwanziger's horizon term, eq.®. This is the task of the present work. 

Before starting the description of the technical aspects, it might be worth to give a short account 
of how the horizon term of the maximal Abelian gauge has been identified. Having constructed 
this term in the Gribov approximation, eq.(@J), our strategy has been that of looking for a pos- 
sible extension of it, which enjoys the properties of localizability and renormalizability. In other 
words, one looks first at a possible nonlocal term which reduces to eq.(HJ) in the quadratic ap- 
proximation. Further, one requires that such nonlocal term can be cast in local form by means 
of the introduction of a suitable additional set of fields. Finally, the resulting local theory has 
to be multiplicatively renormalizable. We underline that the requirements of localizability and 
renormalizability are unavoidable in order to have at our disposal a consistent computational 
framework. These are strong requirements, which yield a very powerful criterion in order to deal 
with nonlocal terms 4 . In the present case, a unique solution for the horizon term in the maximal 

3 In the case of the maximal Abelian gauge, and for the gauge group SU (2), the gap equation defining the Gribov 
parameter 7 reads £Q 



l = l9 



3 , f d 4 k 1 



4 a J (2^) 4 fc 4 + 7 4 ' 

4 See for instance the recent results |39l I4(J| on the construction of a renormalizable nonabelian massive model 
based on the use of the gauge invariant nonlocal operator Tr J d^xF^-^F^. 
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Abelian gauge, which reduces to expression (£Q) and which fulfills the requirements of localizability 
and renormalizability, has in fact emerged from our analysis. It reads 
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where A4 ab stands for the off-diagonal Faddeev-Popov operator of the maximal Abelian gauge 



Notice that the horizon term of eq.flBJ) reduces precisely to expression (jlj) in the quadratic ap- 
proximation, amounting to replace the operator Ai ab by the Laplacian, i.e. by —5 ab d 2 . No other 
solution for the horizon term in the maximal Abelian gauge, fulfilling the requirements of lo- 
calizability and renormalizability, has been found. Although several candidates which reduce to 
expression (j3J) can be easily written down, they are ruled out by the requirement of localizability 
and renormalizabilty. Examples of such terms can be obtained from (j3J) by including more space- 
time derivatives in the numerator, while considering higher powers of the Faddeev-Popov operator 
Ai ab in the denominator 5 . However, the localization procedure of such terms would require the 
introduction of dimensionless auxiliary fields, a feature which jeopardizes the renormalizability of 
the resulting action 6 . Expression (jHJ) can be seen thus as the minimal extension of the nonlocal 
action (£Q). It shares great similarity with Zwanziger's horizon term in the Landau gauge, eq.®, 
as expressed by the presence of the inverse of the Faddeev-Popov operator, (A4~ 1 ) ab , which looks 
rather natural. We remind in fact that the diagonal ghost sector turns out to be unaffected by 
the presence of the Gribov copies, see Appendix B of 0. The proof of the localizability and of 
the multiplicative renormalizability of expression (jHJ) constitutes our main result. 
The present work is organized as follows. In Sect. 2 we give a brief account of the construction 
of the maximal Abelian gauge fixing condition and of the results obtained in PQ on the infrared 
behavior of the gluon and ghost propagators. In Sect. 3 we implement the localization procedure 
for the nonlocal horizon term (jHJ), presenting the BRST invariance as well as the rather rich set 
of additional global symmetries of the resulting local action. Sect. 4 is devoted to a detailed proof 
of the multiplicative renormalizability of the theory. The proof extends to all orders by means of 
the use of the algebraic renormalization [32] • An interesting feature of the theory is represented 
by the existence of a residual local U(l) Ward identity, typical of the maximal Abelian gauge. As 

1 /2 

a consequence, the nonrenormalization theorem, Z g ZJ^ = 1, relating the renormalization factors 
5 An example of a term of this kind is given by the expression J d\, e ah Ma ® Mdc e ch A^. 

6 Even if our present analysis has relied on a direct inspection of possible candidates, it is worth mentioning 
that the requirements of localizability and renormalizability might be linked in a natural way to the deformation 
technique introduced in |41j . which could provide a more systematic way to look at the construction of the horizon 
term for a generic gauge, once the corresponding horizon function in the Gribov approximation has been established. 
Essentially, this would amount to start from the localized form of the horizon term in the Gribov approximation, 
and search for the most general local deformation which preserves all symmetries of the starting model as well 
as its field content. By means of the use of the master equation and of the antifield formalism, all possible 
deformations can be identified in terms of cohomolgy classes of the nilpotent functional operator associated to the 
master equation. The construction of the extension of Gribov's horizon term through the deformation approach is 
under investigation. 




(9) 



D <* = 5 ab d, - ge ab A, . 



(10) 
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of the gauge coupling constant Z g and of the diagonal gluon field Za, still holds in the presence 
of the Gribov horizon. In Sect. 5 we introduce a generalized dimension two gluon operator and we 
prove its multiplicative renormalizability in the presence of the horizon term (JHJ). Furthermore, 
its anomalous dimension is not an independent parameter of the theory, being obtained from 
the renormalization factors of the gauge coupling constant g and of the diagonal anti-ghost field. 
Sect. 6 contains a few concluding remarks on further developments and on potential applications 
for lattice numerical simulations. 

2 The gauge fixing condition for the maximal Abelian 
gauge 

In this section we introduce the gauge fixing condition for the maximal Abelian gauge, providing 
a short account of the main results on the infrared behavior of the gluon and ghost propagators 
obtained in [T|. 

Similarly to the decomposition of the gauge field given in expression (|5jl. for the field strength 
one has 

^ = F^T a + F, u T\ (11) 
with the off-diagonal and diagonal parts given by 

F; v = DfAl-DfA^, (12) 
F, u = d,A u -d u A^ + ge ab A-A b u , 

where the covariant derivative D ab has been defined in eq. (jlUJ) . For the Yang-Mills action in 
Euclidean space one obtains 

S YM = \J d 4 x (f; v F; v + F^) . (13) 

As it is easily checked, the classical action (|13|) is left invariant by the gauge transformations 

8AI = -D; b u b - ge ab A b ^ , 

5A» = -d,u-ge ab A a y . (14) 

The maximal Abelian gauge is obtained by demanding that the off-diagonal components A a ^ of 
the gauge field obey the nonlinear condition 

D*Al = , (15) 

which follows by requiring that the auxiliary functional 

K[A] = J d A xA^Al , (16) 

is stationary with respect to the gauge transformations (|14|). Moreover, as it is apparent from the 
presence of the covariant derivative -Djf, equation (J15|) allows for a residual local U(l) invariance 
corresponding to the diagonal subgroup of SU(2). This additional invariance has to be fixed by 
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means of a suitable gauge condition on the diagonal component A^, which will be chosen to be of 
the Landau type, also adopted in lattice simulations, namely 

= . (17) 

The Faddeev-Popov operator, A4 ab , corresponding to the gauge condition (|15)l is easily derived 
by taking the second variation of the auxiliary functional 7£[A], being given by 

M ab = -D a ^D cb - g 2 e ac e hd AlAl . (18) 

It enjoys the property of being Hermitian and, as pointed out in [IBj, is the difference of two 
positive semidefmite operators given, respectively, by — D a ^D cb and g 2 e ac e bd A c fl A d l . 

As discussed in [T], in order to deal with the existence of Gribov copies which affect the gauge 
condition (|15|). one proceeds by restricting the domain of integration in the Feynman path integral 
to the so called Gribov region C , defined as the set of fields fulfilling the gauge conditions (|T5|). 
(|17j) and for which the Faddeev-Popov operator A4 ab is positive definite, namely 

Co = {A„ A«, d,A, = 0, DfAl = 0, M ab = -D; c D cb - g 2 e ac e bd A c X > 0} . (19) 

The boundary, <9Co, of the region Co, where the first vanishing eigenvalue of M. ab appears, is called 
the first Gribov horizon. The restriction of the domain of integration to this region is supported 
by the possibility of generalizing to the maximal Abelian gauge PQ Gribov's original result ^7j 
stating that for any field located near a horizon there is a gauge copy, close to the same horizon, 
located on the other side of the horizon. 



The restriction to the region Co can be implemented through the no-pole condition on the off- 
diagonal two point ghost function. To the first order, this condition amounts to the introduction 
of the nonlocal term given in expression (0J). As a consequence, the tree level gluon and ghost 
propagators get deeply modified in the infrared. More precisely, both off-diagonal and diagonal 
transverse components of the gluon propagator turn out to be suppressed in the infrared [lj. The 
diagonal component of the gluon propagator is found to display the characteristic Gribov type 
behavior, i.e. 

(Mk)M-k)) = ^7 (V - ^) • (20) 

The off-diagonal propagator turns out to be of the Yukawa type, being given by 

(A«(k)<(-k)) = * a V^(<V-^) , (21) 
a, b = 1,2 

where m denotes the off-diagonal dynamical mass originating from the dimension two gluon con- 
densate (A a Afy As already remarked, the behavior of the transverse diagonal and off- 
diagonal gluon propagators, eqs. fl20|) , (j21j) . is in qualitative agreement with the lattice results [Tfij . 
In the case of the ghost propagator, it turns out that the off-diagonal component exhibits infrared 
enhancement, according to 

S(*)ljfc=o « (22) 

s(*o = ^E^) ca (- fc )) . 
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where (c a ,c a ) stand for the off-diagonal Faddeev- Popov ghosts. Also, the diagonal component of 
the ghost propagator turns out to be not affected by the restriction to the first horizon. 



3 Localization of the horizon function and symmetry con- 
tent 

3.1 Local action from the horizon 

In this section we describe the localization procedure for the horizon term. Let us start by 
considering the partition function of Yang-Mills theory quantized in the maximal Abelian gauge, 
in the presence of the horizon term of eq.flSJ), namely 

Z = J DA a DADb a DbDc a Dc a DcDc exp (— (SyM + 5mAG + SHorizon)) ; (23) 

where Sym stands for the Yang-Mills action, eq.(|I5jL and Smag denotes the gauge fixing term 
corresponding to the gauge conditions of eqs.(|I5 j) .(jr7 j) . The fields (c a , c a ) are the off-diagonal 
ghosts and antighosts, while (c, c) denote the diagonal ghost and antighost. Also, (b a , b) are the 
off-diagonal and diagonal Lagrange multipliers enforcing conditions (fTK|) and (fTTj) . Smag is given 
by the following expression ^Tj 

Smag = J d 4 x (b a DfA b M - c a M ab c b + ge^cD^A^ + bd^ + cd^c + gs ab A«c b )) , (24) 

and ^Horizon ^ s ^ ne horizon term (JHJ) 

^Horizon = 7 V / ^X A, (AT ^ ■ (25) 



As in the case of the Landau gauge fHl ED] , the horizon function (|23|) can be localized by means 



of a pair of complex vector bosonic fields, , <#f ), a, b = 1, 2, according to 



e- 5 — = J D0 J D0(detA^) 8 exp(- J A (^ 6 M oc < + 7 V1< - 4>f) A S) > ( 26 ) 

where the determinant, (det .M) 8 , takes into account the Jacobian arising from the integration over 
the fields (<pff,4ff). This term can also be localized by means of a pair of vector anticommuting 
fields (u^ , oj^) , namely 

(detM) 8 = J DujDuj exp (^J d 4 xufM ac uj^j . (27) 

Therefore, we obtain a local action which reads 

Z = J Dm e ~ SLocal , 
Dm = DA a DADb a DbDc a Dc a DcDcD(j)D(f)Du;Duj , (28) 
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where 

^Local = SyM + 5mAG + Sfa + S 1 , (29) 

with Sfa, S-y given by 

= J d 4 x {^M ac <p cb - ofM ac u;f) , 
S, = J dW* A M . (30) 

3.2 BRST invariance 

In order to establish the BRST invariance of the resulting local theory, we proceed as in [T9"l |2"U] 
and consider first the particular case when 7 = 0, i.e. 

<SLocal|7=0 = ^YM + ^MAG + S<t>u ■ (31) 

In this case we have in fact introduced nothing more than a unity written as 

J D<PD4>DujDuj exp {- J d 4 x {^ ab M ac (j) cb - tf b M ac u cb )^J = 1 . (32) 

Furthermore, the action (}3*Tj) may be written in a BRST invariant fashion. To see this, let us 
introduce the following nilpotent BRST transformations 



sa; 

sA^ 


= -(d^c + gs ab Ay) , 




sc a 
sc 


= ge ab c b c, 
= ^e ab c a c b 
= b a , 






sc a 


sb a = , 




sc 
suf 


= b, 
= < 


sb = , 
sutf = , 

s r b = o , 


(33) 



s 2 = . (34) 

Now, let So be the action defined by 

S = S YM + sj d 4 x {c a D ab A\ + cd,A, + o; b M ac ^) , (35) 

which satisfies 

sS = . (36) 

Acting with the BRST operator s, and recalling the expression of the Faddeev-Popov operator 
M ab , eq.flHJ), we obtain 

S = S YM + S MAG + S^ + j d A x Q ab T ac <f b , (37) 
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with 



T ab = 2gs ac {d,c + gs de Alc e )D* + ge ab d,{d,c + ge cd A c /) 

-g 2 (e ac e bd + e ad e bc )Al{D c *c e + ge ce A e ^c) . (38) 

Expression (|57j) differs from the action (|3*T|) by the presence of the last term. However, following 
[131201, we may transform S\, oca ]\ 1= Q into So by performing the following shift in the variable 



ab 



ab - (M~ 1 ) ac F cd (j)f , (39) 



whose corresponding Jacobian turns out to be field independent. Thus, the following equivalence 
holds, namely 

Dm e~ So = / Dm e - 5 Locail 7 =o _ (40) 



(44) 



Further, let us consider the term 5* 7 given by (|3U|). One can easily check that the term E ab A ll (f> ab 
which can be written as 

e ab A^ ab = e ah s{A^ ab ) - e ab (d,c + ge cd A c /)u ab . (41) 

Once again, according to [P9*l 120] , we can eliminate the last term of eq. (J41j) by means of the change 
of VcLrictblcs 

uf - u ab - (M- l ) ac 1 2 ge cb {d ll c + ge de A d c*) . (42) 
Therefore, for the partition function we obtain the final expression 

Z = J D^-lSo+fs/A^t^tf-fV; 4 ))]. (43) 
Notice also that, due to the identity 

7 2 2 J d*xe ab {A^f ~ s {A^f)) = 7 2 J d'x (D*<ft - a (^ 6 <)) , 
expression (|4*3j) becomes 

Z = J OTe -[*+i 2 /A(«' s (W))] . 
Nevertheless, due to the term D® b <f) b ® , the action 

S + 7 2 J d'x {D^; - s (D; b u b ;)) , (46) 

is not yet BRST invariant, a point which can be dealt with by means of the introduction of a pair 
of BRST doublets of local external sources [111201, (Ufi, M$) and (V^ b , N$), which transform 

as 

suf v = -m;1, sM»l = o, 

*V$ = N% aN* = Q. (47) 



(45) 
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As pointed out in [HI 1201, the introduction of these external sources allows us to promote expres- 
sion pUj) to a BRST invariant action. In fact, let Sources be the action 

Sources = 8 j dk (-U£D*<fc + V£D*U%?) 

d 4 x (m;z Df^ + u;: a (rtftf) + n% D?up + a [Df^) ) , (48) 

which obviously satisfies 

■S*Ssourccs • ( / 1^) 

Moreover, when the sources U, M, V, N attain their physical value [TH| I20|. defined by 

1V1 fiv\phys — V ^ v Iphys — 0^7 , 

k^Jphys = N fu Iphys = , (50) 

it immediately follows that 

Ssourceslphys = J tfltfft ~ >? 8 («)) , (51) 



which is nothing but expression ([46)1 . One sees thus that the use of the external sources U, M, 
V, N enables us to introduce an extended action E 

^0 = So + S'sources , (52) 

which enjoys the important property of being BRST invariant, 

sE = , (53) 

while reducing to expression when the sources attain their physical value, eq.(|50jl . Explicitly, 
we have 



Sym + S MAG + s J dh «AT C < ~ U£ Dftf? + V£ Dfu b :) 

S YM + S MAG + J A ( r b M ac rf - ufM ac u;£ + tf*F*<ft + M£ Df^ + N£ Dfu, 



+U™[D a y v c + 9e a \d,c + ge^A^)^} + V™[D a ffi + ge ab (d,c + ge de A*c e )u b u c }) . 

(54) 

3.3 Inclusion of the quartic ghost term 

Although being BRST invariant, the action S is not yet the most general classical action to 
start with. As discussed in previous works (TTJ H2] , the nonlinearity of the gauge condition (JT5J) 
requires the introduction of a quartic term in the Faddeev-Popov ghost fields 

VcVcV , (55) 
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which is in fact needed for renormalizability purposes. The parameter a in expression (J55j) is a 
gauge parameter. In our case, due to the presence of the localizing fields (0^f, <^f, cu^ b , uj^ b ), the 
quartic ghost term ()55|) is introduced in a BRST invariant way through the following action E a 

E a = s j d A x^ {c a b a - ge ab ^c + g 2 ^ (fiftf ~ - 2g 2 ^^c b c b ) 

J d A x (b a b a - 2ge ab b a c b c + g 2 ^c a c b c b + g 2 ~ «) ~ ^M d ) 

-2g 2 {4%4$ - u a ;u a ;) J>c b + 2g 2 ^^b b c b - 2g^<j>fe bc ^c c c) . (56) 



a 
2 



Adding the BRST invariant term E a to the action E , modifies the equation of motion of the 
off-diagonal Lagrange multiplier b a , according to 

* (S °^ Sa) = D tK + « ( ba ~ 9e ab c b c + g 2 ^c«) . (57) 

Therefore, we see that maximal Abelian gauge condition (|15p. D^A b = 0, is attained in the limit 
a — ► 0, which has to be performed after the removal of the ultraviolet divergences JHE]- We 
also remark that expression (|5fijl contains a unique free parameter, namely, the gauge parameter 
a. As a consequence, the whole term E a vanishes in the limit a — > 0, allowing us to integrate out 
the auxiliary fields (<^f , <^f , tu® b , u)^ b ), thus recovering the horizon term (JSjl. As we shall see in the 
next section, this important feature follows from the fact that the action (So + E a ) fulfills several 
Ward identities which, in particular, uniquely fix the form of the term E a . 



3.4 The global U(8) symmetry 

In addition to the BRST invariance, and in complete analogy with the case of the Landau gauge 
[THU^D], the action (E + E a ) displays a global symmetry U(f), f — 8, expressed by 



Qt (So + S Q ) = , (5* 



with 



■I'" I ' \ l' A, ,<■'< " Jf^ca 1 M facb v fijjjca 

r c c X \ 

_i_T/ ca ]U~ cb i AT ca _ TT cb | (^Q) 

alM SV cb ° u 5M ca a ^5N cb au SU ca ' 

The presence of the global invariance U(8) means that one can make use [TTHI27I] of the composite 
index i = (a, /i), i = 1, . . . , 8. Therefore, setting 

«^?,«) = (0?,^,wf,cDf), (60) 

and 

Kl, v; b , m«1 Np = {v%, v%, m; 4 , at;j , (ei) 
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A 


b 


c 


c 





(j) LO 


uJ 


u 


V 


M 


N 


dimension 1 


2 


2 





1 


l l 


1 


2 


2 


2 


2 


gh number 





-1 


1 





1 


-1 


-1 








1 


<2/-charge 











1 


-1 1 


-1 


-1 


1 


-1 


1 



Table 1: Quantum numbers of the fields and sources 



we rewrite expression (p)2"|) as 

S + S Q = Sym + W + J d'x {&M ab tf - ^M ab u b + Q*F*<ft + M;. Df<j>\ + iV» 

+U^[D ab u b + ge a \d,c + ge^A^tf] + V^D^ + gs ab (d,c + ge ca %c d )cu b ] ) 
+ | J d 4 x (b a b a - 2ge ab b a c b c + g 2 c a c a c b c b + g 2 - u?u?) (4>)<j>) - u b u b ) 

-2g 2 - Qfuf) c b c b + 2g 2 uj*<j) a i b h c b - 2g 3 ^^e bc c b c c c) (62) 

For the symmetry generator we have 

5 T „ S S S 



Qii = d% 6?—-tf-=-+L0?—-L0 



A A A A \ 

By means of the trace of the operator Qy, ie., = Q/, the z-valued fields turn out to possess 
an additional quantum number, displayed in Table ((H), together with the dimension and the ghost 
number. 



4 Identification of the final complete classical action E 

Let us proceed by establishing the rich set of Ward identities which will enable us to analyse the 
renormalizability of the theory to all orders. Let us first identify the final complete action to 
start with. To this purpose, we have to properly define the composite field operators entering the 
BRST transformations (j33J),((I7J). We notice that the BRST transformation of the gauge field A® 
can be written as the sum of two composite operators, i.e. 

sA; = O x + 2 (64) 

where 

1 = -D; b c b , 2 = -ge ab A\c . (65) 
Moreover, thanks to the nilpotency of the BRST operator, s 2 = 0, it follows that 

sOx = -s0 2 ■ (66) 

Therefore, the two composite operators, 01,02, can be defined by means of the introduction of 
the external sources (0",t",£"J 

E« = j d*x (n; {-Dfc b ) + r; (-ge"»A»c) + £ s(-ge ab A,c)) . (67) 
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To guarantee the BRST invariance of S^t, we require that 

s^l = st; = 0. (68) 
Moreover, it is easily checked that the new action 

S + S a + E« , (69) 
is left invariant by the following set of transformations: 
• the 5i symmetry 



the Si symmetry 



S t c a 


= €, 


5$ 


= S-c a 


S t b a 


= ge*<fic 


mi 


= V£, 


~S l c a 










= ge ab u b c 




= -U a - 



(70) 



(71) 



with 



Si (e + S Q + £«) = ~S t (S + S Q + E« ) = . (72) 

As transformations (f70|) . f7T{) contain composite field operators, i.e. ge ab (f) b c and ge ab uj b c, we 
define them by means of external sources (77", A") and giving rise to two set of BRST 

doublets 



so that 



< = A", S A° = 
^ = pi sp« = 0, (73) 



= y (ge ab \y b c + 77? [^£ afe ^c + 9 - e ah e c %\c c c d \ 

+ge ab p a i uj b c - ^[ge ab $c - 9 ^e ab e cd u b c c c d ] ^ . (74) 
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Further, introducing the BRST invariant sources (f2 M , L a , L), coupled respectively to the nonlinear 
BRST transformations of the fields A^, c a , c, for the final expression of the complete action E we 
shall start with, we have 

E = S + E Q + E« + E<2 + s J d 4 x {-%A, + L a c a + Lc- X U%V$ , (75) 

namely 

X = S YM + Smag + s J d A x (QfM^t ~ Ki D % + ^ ^ + ^"V^c + ^i??wf c 
-n^Al - ge a %A\c - %A, + L a c a + Lc- X U*y$ 
+s J <£ X °L (Tb a - ge ab ^c + g 2 ^ - u> b u b ) - 2g 2 ^^c b ) 

= S YM + J d A x (b a DfAl - c a M ab c b + ge ab c a cD b ;Al + bd^ + cd^c + ge ab A a /) 

+^M ab <p\ - u%M*ui + fl^tf + M« Dftf + UftDfit + ge ab (d»c + ge cd A c /)^ 
D ab u b + V£[D*ft + ge ab (d,c + ge cd A c fl c d )oj b } - Q; Dfc b - ge ab r^c 

2 

+Q[ge ab (D b ll c c c )c - 9 —e ab e cd A\c c c d \ - Q^c + ge ab A a /) + ge ab L a c b c + 9 -e ab Lc a c b 

2 2 

+ge ab \ a i( j) b c + r,*[gs ab u b c + 9 je ab e cd (j) b c c c d } + ge ab p a t Gj b c - ^[ge ab $c - 9 —e ab e cd Q b c c c d \ 

+ x(AOT + djjvj) ) 

+| J d 4 x (b a b a - 2ge ab b a c b c + g 2 c a c a c b c b + g 2 - Qfuf) - u b cu b ) 

- 2g 2 - ufuf) c b c b + 2g 2 uo a l (j) a i b b c b - 2g 3 u^s bc c*c c c) . (76) 

One should notice that the term quadratic in the sources (M" i; V^, U^, Nfy) , i.e. x(M^V^ + 
UuiN^), in eq.(|75jl is of dimension four. As such, it is allowed by power counting and has to be 
added for renormalizability purposes. The parameter x stands for a free coefficient. 

Expression (|7B|) represents our starting action. 



4.1 Ward identities 

It turns out that the classical action S, eq.(!76|). obeys the following set of Ward identities: 
• the Slavnov- Taylor identity 



5(E) = , (77) 



with 



I ,] i i 5S 5E\ 5E 5E 5S 5E 5E 5E5E , „ 5E ,6E 
5(E) = / cfx — + ———+ + — — + — — + b a — + b— 

1 ' 1 <Jft£ St" J 5A a Skip 5A^ 5L a <5c a 5L 5c <5c a 5c 



5E , xa 5E , ^ a 5E\ 



(7f 
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the Ward identities corresponding to the Si and Si symmetries, eqs.flTOj) , (J7TJ), i.e. 

W((E) = J d'x (« g + , g + vsg - * § + g§) = , (79) 

Jv. (E ) S / A (* g - c'g - ^ + , ? § + g§) = , (80) 

the Qj(E) and Qi(S) identities obtained by commuting the Slavnov- Taylor identity, eq. (|77jl , 
with the Wi(E) and VVj(£) identities, eos. (|75)> . (IKU1) . namely 

Qi(S)=0, (81) 



, v-. / ,. / / ^ n \ ^ ST f ST a \ ST ST irn ST _ ST 

Qi(T) = / — + u$ — + — — + - 77^ + ^ - C77^ + Ki' 



SXf 1 ) Sc a S^Sb a \8$ ri JSL a Suf ^Sn a 



and 



Qi(T) = , (82) 

a(E) s , a 1 1 _ _ «j _ + __ + ^_ _ ^ j ^ - ^ + ^ - m» - 

the rigid 72.— identities 



*g>(E) S / A (tf g - ^ + ^ + ^ + *g - tfg) = , (8 



3) 



K t] (T) = y ^ _ - + ^ — - ^ _ + - M M - 



+ 7V a 

1 * ///, 



"» SN a - 



(84) 

«« m - / A (* g - *g - ^ - I* * - *g + *g) - , M 
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The trace of Tlfj, eq.(JE3J, defines the Q/-charge of all fields and sources. 
• the diagonal gauge-fixing condition and the antighost equation 



5S 
~8b 



d»A„ , (87) 



8E _ 5E 



8c "sn 



the £>(E)-Ward identity 



the local C/ (1) Ward identity 



with 



<5c a 5L a 



W 3 (S) = -d 2 6 , (90) 



p r c r A A A \ 

+M a h 7V a h77 a h$ a h A a ho" h L a I (91) 

+ *» 8M% + V 8N b t + Vl 8 V b t + 1 5$ 1 8X\ + Pl 8p\ 8L b J ' [ * } 

Notice that the breaking term in the right hand side of eq. (l9U|) . i.e. d 2 b, is linear in the 
quantum fields. It is thus a classical breaking, not affected by the quantum corrections 

4.2 Useful commutation and anti-commutation relations 

For further use, let us give here the relevant commutation and anti-commutation relations between 
the nilpotent linearized Slavnov- Taylor operator B^, 

£s£s = , (92) 
.,.,,££ <5 8E ( 8 5 \ 8E 5 8E 8 



8Q* 5r«y 5A« <L4£ \6n* 5t«J Stt^SA^ 8A fl 8tt IM 
8Z 8 8Z 8 8 8Z 8 ,J ,5 5 

+ fr - „ + - fr + 777- + -7-77- + 6 a 7— + 6— + t— 

5L a 5c a <5c a <5L a 5c 5c 5L 8c a 8c 8$ 

+*i + ~ M -sk ~ ^ ~ T » ] h + K i + O • (93) 
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and the linearized operators Wf, VVf , Qf, Qf, Ps corresponding to the Ward identities (|75j). 
PH. flEU), (E2D, HP, given by 

VK J X \ W * 5c- C Suf ^ 5Q- + Vi 5L° + 5p? 5b- + 5b- 5p? J ' 
^e [ * (fSZ a \ 5 5Y, 5 5T, 5 5T, 5 ( 5E 

Qf = (fx — + -7— ft" + , „ + "Ti — ? h ( v Pi 



^ S + ^-kA + K^) , (96) 



Qf = d 4 x —-$)— + — — + — — + — — + — - XI 



\5pf ) 5c- 5c a 5p a i 5^5b a 5b a 5ti? \5uf 1 J 5L a 



A _ e — + U a -—- M a — ) (97) 



and 



£> s = d A x [ c a — + Yi~Tr~ + Ti — ft - • 98 



The commutation and anti-commutation relations of these linearized operators are found to be 

{VVf, Be} = Qf, 
[Wf,S s ] = Qf, 

[Pe,B e ] = 0. (99) 



Furthermore, we also have 



r 5 „ , 5 



[W 3 ,£ s ] = 0. (100) 



5 Algebraic characterization of the most general counter- 
term 

According to the algebraic renormalization [12], the most general local counterterm Sqt allowed 
by the Ward identities (JZZD , CSD, (IHD, (ED, ©, (El, ©, (ED, (JHZD , (JHHD , (H, O, is 
an integrated local polynomial in the fields and sources of dimension four, see Tables (J2J) and (|5jL 
which obeys the following set of constraints 

i3 E S CT = 0, (101) 
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A 


b 


c 


c 





(j) UJ UJ 


dimension 


1 


2 


2 





1 


1 1 1 


gh number 








-1 


1 





1 -1 


Q /-charge 














1 


-1 1 -1 



Table 2: Quantum numbers of the fields 





n 


T 


t 


L 


M 


N 


U 




7? 


A 





P 


dimension 


3 


3 


3 


4 


2 


2 


2 


2 


3 


3 


3 


3 


gh number 


-1 


-1 


-2 


-2 





1 


-1 





-2 


-1 


-1 





Q/-charge 














-1 


1 


-1 


1 


-1 


-1 


1 


1 



Table 3: Quantum numbers of the sources 



(102) 
(103) 



(104) 

(105) 
(106) 

0, 

. (107) 

From eqs. (H07|) it follows that Sct is independent of the diagonal Lagrange multiplier b, and that 
the diagonal antighost c enters only through the combination (f2 M + d^c). Furthermore, from 
general results on the cohomology of gauge theories EH] j it turns out that the most general 
solution of the constraint (jlOlj) can be written as 

Sct = a S YM + B S A^ , (108) 
with A*^ 1 ) being an integrated local polynomial with ghost number -1, given by 



WfEcT = 0, 

HfS CT = 0, 

QFSct = o, 

Qf Sct = 0, 





= 




= 




= 


Sct 


= 



T^t, Sct — , 
W 3 S CT = , 



8b 

^Sct n. <^Sct 
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A ( - 1} = J (fx {a x Q;A; + a 2 r^A" + a 3 Q ge ah A\c + a 4 Q d,c a + a 5 ^ ge ab A,c b + a 6 (d^A* 

+a 7 (n^ + fyc)^ + a 8 c a L a + a g cL + a 10 77? ge^c + a n rjfuf + a 12 tf a ge ab u b c 
+a 13 + a 14 + ai 5 rfa>? + a 16 E/£ + a 17 L£ + a 18 V;, d^u? 

+a 19 V£ 0£%utf + a 20 c a b a + a 21 0e°Wc + a 22 ge^A^ + a 23 0?$$$ 
+a 24 cDf #u>X + a 25 ^^05 + a 26 cDf <^ a a;j + a 27 + a 28 uttftfrf 

+a 29 + a 30 + a 3 i wf^w} + a 32 cD^c 6 + a 33 utf^cV 

+a 34 cDf + a 35 flftfA^ + a 36 + a 37 ^<^4 + «38 

+a 39 w? ge^AMi + a 40 cD« ^(c^)^ + a 41 x^^) , 

(109) 

where the coefficients a*, i = 0, ...,41 are free dimensionless parameters. Notice also that in the 
derivation of expression ()1U9|) use has been made of the fact that the action £, and thus £ct, are 
left invariant by the following discrete symmetry 

#i _> #i j 

^2 _^ _^2 ^ 
^ydiag _\J/diag ^ /-qq\ 

where \l/ a , a = 1, 2, and \l/ diag stand, respectively, for all off-diagonal and diagonal fields and exter- 
nal sources. As one can easily recognize, this symmetry plays the role of the charge conjugation 

nu. 

After a rather long calculation, from the constraints ([102]) . ([103]) . ([104]) . ([105]) . ([106]) . one finds 
that 

dy = (2io = a ll = a 12 = a 13 = = 0-15 = , 

&25 = 0-26 — a 27 — 0-28 = a 29 = a 30 = a 31 = > 

^33 = ^34=0, (111) 

while 



«3 
«5 

ai6 

«23 
«32 

a 3 s 



ai , 
-a 4 , 
a 6 , 

— CLyj — — fl^g = CZ 4 g = CL\ — fl 4 -f- 
—(120 j 

-a 24 = # 2 a 20 - a# 2 a 8 , 
-2g 2 a 20 + ag 2 a 8 , 



-a 4 o 



^39 

~2 



a 35 a 3 6 



(137 



r 



(112) 
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Making use of the relation 

t ;a; = q;a; + ftrtfj + g^Afr) + b^a;) , (113) 

and remembering that £>| = 0, we obtain 



A (_1) = 



J d'x ( + a 2 )plAl + ^^c] + (a 2 + a^D ab c b 

_ a 4 + a 8 )[^V- - ^<^ a ] - K + as) c^ATV- 
-a 6 c a Df A\ + a 8 L a c a + a 9 Lc + a 41 xU^V* + a 20 (c^ - ge ab c a c b c) 
- (aa 8 - a 20 ) g 2 ^ - tfrf) + («a 8 - 2a 20 ) g^WtfJ) (114) 

Finally, by renaming the coefficients as follows 

ai + a 2 — > ai , 

a 2 + a 4 — > -a 2 , 

&8 ~~ * 03 , 

^6 — > a 4 , 

^9 — * «5 > 

a 4 i — ► a6 , 
a 

a 2 o -> 7^7 (115) 
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for the most general local allowed counterterm £ct, we obtain 

'1 



^CT 



-Xd,A a u )(d,Al - d u AD - ge a \d,Al){A,A b v - A V A 



I 6 ) 

i-i> 



9 



+ge a (d^A^A^A^ + —(A^A^A^A® + A^A^A^A^) 
+(a + Aa^AlAlAlAl + (a, - a^D^A* 



a 



^^(d^A^Av-dvAJ 



-(a 3 + a 4 )[c a dV - cV^aAK - 2eV% V - aVcM^] 
+ (a 3 + a 4 )[^a 2 0^ - tfge a \d,A,)4>\ - 2tfge ab A,d,4> b - g 2 4>^A,A,} 
-(a 3 + a,)ptd 2 ut - u?ge ab (d^)u b - 2^ge ab A^u b - g 2 ^A,A,\ 



+ (2ai - o 3 - a 4 )g 2 e ac e bd {c a c b - ^\ + wfw^Aj 

_a?>— a „ t-)6c a c 



-(a a - 04 - a 5 )0e < VcD*A° - ( ai - 03) + ^cJ^A^ 
-(ax - 2a 3 - a A )[2g 2 e ab e cd u a l Alc d d l 



adj}c\ 



+ g 2 e ab s cd Otd,(Alc d ) ( j ) : 



+2g*e bc u?tfA ft A b l c c - a '(e<V d + £ a V>^(£-c e )$] 
+ (2ai — a 3 — a 4 — as)o 3 (5 



ae^.bd 



5 be e ad ) cotA'A^c 



-(aa + 2a 3 )s + V^})^ + (a 4 + a 2 + a 3 )fi°Z^c b 

-a 2 r;D^c fc + a 5fl e a V;4c - (a 4 + a 2 + a 3 + o 6 )ye a6 ^(D?c c )c 



, ^ab^cdca \b n c n d „ n ^abjaJ>„ ln„ „ \@ ^0,1 ' J^^fl ^ 



+ (a 2 + 2a 3 )^£ aft £ C£i C<c c c d - a 5 2£ a6 LVc - (2a 3 - a 5 ) V- 



+a 5 (fy + fyc^c - a 3 a + ^)c c c d - a^ge ab {\ a i 4> b + ^ 

0?#)c - (a 3 + a 4 + a 5 )[2^f (d,c)d^ b + ge^d 2 ^ 



+P^t 



-■'AA,<\ 



{ ai + a 2 + a 3 + a^ge^^m 



flirt 



fit I , 



+2g z u, 

_( ai + a2 + a 3 )(M;^«V- + Zftitfw? + V^D ab (t) b + ^2#0f) 



13 ) c a c a c b c b 



-a 6 x(M«v; i + [/«Ag) 

+ | y d 4 x (a 7 fo a 6 a - 2 (07 - a 5 ) ge ab b a c b c + g 2 (a 7 - 2a 

+ a 2 (07 - 2a 3 ) (^ - - *>M) " 2 ^ ^ " 2a 3) - ©X) ^ 

+ 2o 2 (a 7 - 2a 3 ) w?^6 6 c 6 - 2o 3 (a 7 - 2a 3 - a 5 ) ^°£ 6c c 6 c c c) . (116) 



5.1 Stability of the classical action and renormalization factors 

After the characterization of the most general local counterterm £ct ; eq. ()116jh compatible with 
all Ward identities, we still have to check that it can be reabsorbed through a multiplicative 
redefinition of the fields, sources and parameters of the starting action E, according to 



where 



S(^o,^o,^o,fio,r ,Co) = S(^,^,J,fi,r,C) + eE CT (^,^,J,fi,r,C) + 0(e 2 



^0 = 4 /2 * , 



(117) 
(118) 
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denotes the off-diagonal fields, \I/ = (^4°, b a ,c a , c a ), while 



^ = 4 /2 ^' ( 119 ) 

stands for the fields, i/j = [A^, b, c, c, </>", 0", , u) 4 a ) . 
Also J and ( 

Jo = Zj J , 

Co = ^ f C, (120) 

denote the external sources, J = (£" L a , L, f^, J7"j, V™ , M^, iV^, 77?, A", pf) , and the parame- 
ters ( = (g,x,ct), respectively. 

Moreover, recalling that the sources and are coupled to composite operators, 0\ = —D^d 3 , 
2 = —ge ab A b c, displaying the same quantum numbers, see eqs.([65 |) .()67j ) . they renormalize as 

C9 = ** • (i2i) 

where the matrix 



allows for the mixing at the quantum level between the operators 0\, 02. By direct inspection of 
Sct, the renormalization factors are easily found to be 



z b 




z]z- c 'z c , 


z b 




7 2 
Z 9 ' 


Zn 










Zg 2 ZcZ c , 


Zu) 




z 2 a z- c x z c , 


Z<f> 




z§ = z c , 


Zu 




7— 1 /71/2 7 

Z g Z c Z ^ > 


Zn 




7 (7— 1/2 

Z/ 9 Z/ g Z/y , 


Zm 






Zf\ 




Z g 2 Z- c Z^ 2 , 


Z\ 




7 7— 1 7I/2 


z P 




£-1/2 








Zl°- 




z- l zl' 2 z;^ , 


Z L 




Z -l Z l/2 Z -l/2 ; 


Za 




7 -2 



(123) 
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with 



Finally 



Z A = 


1 + e (a + 2ai) , 


z 9 = 








Zc = 


Z 5 = 1 - e (a 3 + a 4 ) , 


z c = 


1 + e (a 3 — a 4 — 2a 5 ) , 


— 


1 c In n \ 


Zy = 


1 + e ^— ai - a 2 - — + 


z x = 


1 + e (2ai + 2a 2 + 03 — 04 


Z a = 


1 + e (07 + ao + 2a 4 ) . 




/7 7-1 7I/2 7 -i 7-1/2 
/ Zjy Zj g Zj z Ziy l a 


Zq t = j 






{ z-^z^z-^ 



a4' 

T 

04 — ae) , 

(124) 



(125) 



Equations ()123|) , ()124j) , ()125J) show that the counterterm S CT can be reabsorbed by means of a 
redefinition of the fields, sources and parameters of the starting action £, establishing thus the 
multiplicative renormalizability of the theory. Let us end this section by noting that the renormal- 
ization factors (Z^, Z^, Z^, Zj) of the auxiliary localizing fields (</>", 4>f, uf, uf) are not independent 
quantities, being expressed in terms of the renormalization factors of the gauge coupling constant 
and of the Faddeev-Popov ghosts. Again, this feature stems from the Ward identities ([79]). (jSUj) . 

We also notice that the renormalization factor Z g of the gauge coupling constant g can be ob- 
tained directly from the renormalization factor Za of the diagonal component of the gauge field 
Afj,. Thus, the nonrenormalization theorem of the maximal Abelian gauge [TH IT2*] 

Z A Z] = 1 , (126) 

remains valid in the presence of the Gribov horizon. This important property is a direct conse- 
quence of the local U(l) Ward identity ([9(J|h 

6 Introduction of a generalized dimension two local oper- 
ator 

As remarked in the Introduction, the inclusion of the horizon term, eq.flSJ), does not prevent us 
from defining a local composite dimension two operator 0^42 

A * = \AIAI + a (c a c a - m + ufuf) , (127) 

which generalizes the gluon operator (^A^A® + ac a c a ) already considered in the maximal Abelian 
gauge [H] , and proven to be renormalizable to all orders ^2] • From the equation of motion of the 
Lagrange multiplier b a 

^ = DfAl + a(b a - ge ab c b c + g 2 ^c a ) , (128) 
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it follows that the integrated operator j d 4 xO 'a 2 enjoys the property of being BRST invariant 
on-shell, namely 

s j d i xO A 2 = j d A x (- {D ab c b ) Al + ab a c a - agc a e ab c h c) = j d 4 xc a ^ . (129) 

As proven in the Landau gauge [13J EHJ 1221 m the case of the operator A^A^, A = 1, N 2 — 1, 
eq. ()129|) implies that the local composite operator Oa 2 is multiplicatively renormalizable. This 
can be established by introducing Oa 2 in the starting action S, and by repeating the algebraic 
analysis done in the previous section. More precisely, we introduce the operator Oa 2 by means of 
a BRST doublet of external sources (J, A) 

sA = J, sJ = 0, (130) 

so that 



£ = £ + s j d A x (AO A 2 + |A J 

= £ + J d 4 x ( JO a 2 + A (A a ^Dfc b - ab a c a + ag^e ab c b c) + | J 2 ) . (131) 

The parameter g is a free parameter, needed in order to account for the ultraviolet divergences of 
the correlation function ((A^(x)A^(x)) (A b v (y)A b l/ {y))) 



The extended action £ is easily seen to obey the following Slavnov- Taylor identity 

<S(£) = , (132) 

with 



61 ^ ' ~ ' ' ' 1 1 5Q° + 6^) 6A* + 60^ III + Jh a 5c? + II Id 

, a SE ,<5£ a - Arfl „ rn 

+6 a — + b— + uf— + 6f— + N°;- —— - M" 



5c a 5c *<50f ^<5u)f "W^ ^SU^ 

" T ^ + A ^ + ^ + J 5A ) ■ (133) 

All other Ward identities, eqs. f|79l — 190 |) . turn out to hold, remaining unmodified by the introduc- 
tion of the sources ( J, A) . Moreover, there is an additional Ward identity 

tf ( v) = / ^ I gl + c^ I =0 . (131) 



5A Sb a 

expressing the fact that the integrated operator j d 4 xC>A 2 is BRST invariant on-shell. 

By repeating the same analysis as before, for the most general local counterterm we obtain 

S CT = a S YM + #g A^ 1 ) , (135) 
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with 



A ( - 1} = J A (ai [Q;A; + C»ge ab Alc} - a 2 $D?<* + a 3 L a c a - a 4 ^D^A\ 



+a 5 Lc + a 6X UX + (a, + a 2 + a 3 ) [MW - V^Dfu\] - (a 3 + a 4 ) utfA^? 



+ -a 7 (c a 6 a - ge ab c a c b c) + - (a 7 - 2a 3 ) <7 2 ^ a - c^ 6 ) - a (a 7 - a 3 ) faWtfJ 
+A Q(og - a 4 ) + a (07 + a 3 ) c a c a - aa 7 - u^) + a 8 | J H . (136) 

In particular, for the renormalization of the sources J and A, and of the parameter g one finds 



— Zj J — ZgZ 5 1 J , 

-1/2 



A = A = Z g Z^ A , 
£> = -Z'e Q = (1 + e(2a - 2a 3 + 2a 4 + a 8 ) ) g . (137) 

We see thus that the renormalization of the source J, and thus of the composite operator coupled 
to it, i.e. O^z, can be expressed in terms of the renormalization factors of the gauge coupling 
constant g and of the diagonal antighost c, meaning that the anomalous dimension of Oa 2 is not 
an independent parameter of the theory. Again, this result is in complete analogy with the case 
of the Landau gauge, where the corresponding operator A^A^ is multiplicatively renormalizable 7 
[22] in the presence of Zwanziger's horizon function (JHJ). 



7 Conclusion 

In this work we have pursued the study of the maximal Abelian gauge by taking into account 
the restriction of the domain of integration in the Faddeev-Popov quantization formula to the 
Gribov region Co- Such a restriction is needed due to the existence of Gribov copies. The Gribov 
approximation, previously introduced in pQ, has been improved through the introduction of the 
nonlocal horizon function, eq.(JBJ). As in the case of the Landau gauge fOJESI!, the horizon term 
of the maximal Abelian gauge can be cast in local form with the help of additional auxiliary 
fields. The resulting local action S, eq. ([76)1 . turns out to be multiplicatively renormalizable to all 
orders. This is the main result of the present article. It could open new perspectives, motivating 
further analytic studies of the maximal Abelian gauge which might provide a better comparison 
with lattice numerical results fHl IS] ■ A partial list of the topics worth to be analysed can be 
summarized as follows: 

• The possibility of having at our disposal a local and renormalizable action incorporating the 
effects of the Gribov horizon enables us to work out the higher order quantum corrections 
affecting both the diagonal and off-diagonal gluon propagator. We expect that the incorpo- 
ration of these quantum corrections will provide a more direct and reliable comparison with 
the lattice data. As an example, we mention the possibility of investigating by analytical 
methods the behavior of the longitudinal off-diagonal component of the gluon propagator, 
which we were unable to discuss within Gribov 's quadratic approximation Lattice re- 
sults [in] have pointed out that this component is nonvanishing. We remark, however, that, 

7 Here too, the anomalous dimension of A^A^ is not an independent parameter of the theory |22|. 
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from the analytic point of view, it vanishes at the tree level. Nevertheless, it might arise at 
the quantum level due to the contribution of the diagonal and off-diagonal transverse com- 
ponents of the gluon propagator in the loop integrals of higher order Feynman diagrams. 
As these transverse components carry nonperturbative information, embodied in the Gribov 
parameter 7 and in the dynamical gluon mass m, a nonvanishing longitudinal component 
might show up at the quantum level. Certainly, the renormalizability of the starting action 
S is a crucial ingredient here. 

• A second issue to be investigated is the infrared behavior of the off-diagonal ghost propaga- 
tor (c a (/c)c b (— &)), which can be decomposed into symmetric and antisymmetric part with 
respect to the off-diagonal indices a,b — 1,2. As already mentioned, the symmetric part 
of the ghost propagator, i.e. Yl a b^ ab (c a (fc)c b (— has been found to be enhanced in the 
infrared region, in the Gribov approximation pQ, as expressed by eq. fl22|) . A more detailed 
study of the behavior of the symmetric component of the ghost propagator would provide 
an important check in order to establish whether this infrared enhancement will be left un- 
modified by the higher order corrections. 

Also, the study of the antisymmetric part of the ghost propagator, i.e. e ab (c a (k)c b (— k)), 
would provide information about the phenomenon of the ghost condensation in the presence 
of the Gribov horizon. The possibility of the formation of a ghost condensate (e ab c a (x)c b (x)} 
in the maximal Abelian gauge, proposed in [37|, has been investigated by several authors 
[1H1 EH1 EDI ■ This phenomenon has been shown to occur also in the Landau gauge [HUEHEB], 
where a few lattice results are available |5H I55j. Here, the ghost condensation has been es- 
tablished by constructing the effective potential for the operator e ab c a (x)c b (x). The existence 
of a nontrivial minimum for the effective potential resulted in a nonvanishing condensate 
(s ab c a c b ) jSU E2 EH] • It could be interesting to see how the phenomenon of the ghost con- 
densation would be modified by the presence of the Gribov horizon, a feature which can be 
now faced. This might result in a better comparison with the available lattice results 8 on 
the ghost condensation, which are done by considering always gauge configurations which lie 
within the Gribov region. The construction of the effective potential for the ghost operator 
e ab c a (x)c b (x) in the presence of the Gribov horizon is under investigation. Both Landau and 
maximal Abelian gauge will be covered. 

• The multiplicative renormalizability of the local composite operator O^, eq. (|127|) . opens 
the possibility to perform a study of the dimension two gluon condensate (Oa 2 ) when the 
restriction to the Gribov region C is taken into account. This might improve our under- 
standing of the dynamical mass generation for off-diagonal gluons, a feature relevant for the 
dual superconductivity picture of color confinement. 

• Finally, we would like to call attention to the nonrenormalization theorem 

Z A Z 2 g = 1 , (138) 

which still holds in the presence of the horizon. This result, stemming from the local U(l) 
Ward identity could motivate further studies of the maximal Abelian gauge from both 

8 We are indebted to A. Cucchieri and T. Mendes for providing us their preliminary results about the numerical 
studies on the lattice of the ghost condensation in the maximal Abelian gauge. 
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Schwinger-Dyson approach and lattice numerical simulations. Concerning the potential 
applications to the Schwinger-Dyson equations, we remark that eq. fll38|) could be employed 
in order to obtain a useful truncation scheme. Also, it would be interesting to look at 
relation (|138|) from the lattice point of view, as it could allow to study the infrared behavior 
of the running coupling constant by analysing the form factor of the pure diagonal gluon 
propagator. 
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